Metabolic alterations and mitochondrial dysfunction in human airway BEAS-2B cells exposed to vanadium pentoxide.

Vanadium pentoxide (V+5) is a hazardous material that has drawn considerable attention due to its wide use in industrial sectors and increased release into environment from human activities. It poses potential adverse effects on animals and human health, with pronounced impact on lung physiology and functions. In this study, we investigated the metabolic response of human bronchial epithelial BEAS-2B cells to low-level V+5 exposure (0.01, 0.1, and 1 ppm) using liquid chromatography-high resolution mass spectrometry (LC-HRMS). Exposure to V+5 caused extensive changes to cellular metabolism in BEAS-2B cells, including TCA cycle, glycolysis, fatty acids, amino acids, amino sugars, nucleotide sugar, sialic acid, vitamin D3, and drug metabolism, without causing cell death. Altered mitochondrial structure and function were observed with as low as 0.01 ppm (0.2 µM) V+5 exposure. In addition, decreased level of E-cadherin, the prototypical epithelial marker of epithelial-mesenchymal transition (EMT), was observed following V+5 treatment, supporting potential toxicity of V+5 at low levels. Taken together, the present study shows that V+5 has adverse effects on mitochondria and the metabolome which may result in EMT activation in the absence of cell death. Furthermore, results suggest that high-resolution metabolomics could serve as a powerful tool to investigate metal toxicity at levels which do not cause cell death.

Floral resource partitioning of coexisting bumble bees: Distinguishing species-, colony-, and individual-level effects.

Resource partitioning is considered a key factor in alleviating competitive interactions, enabling coexistence among consumer species. However, most studies have focused on resource partitioning between species, ignoring the potentially critical role of intraspecific variation in resource use. We investigated floral resource partitioning across species, colonies, and individuals in a species-rich bumblebee community in the diversification center of bumblebees. We used a total of 10,598 bumblebees belonging to 13 species across 5 years in the Hengduan Mountains of southwest China. First, we evaluated the influence of a comprehensive set of floral traits, including both those related to attractiveness (flower color and shape) and rewards (pollen, sugar ratio, nectar volume, sugar concentration, and amino acid content) on resource partitioning at the species level in bumblebee-plant networks. Then, we explored intraspecific resource partitioning on the colony and individual levels. Our results suggest that bumblebee species differ substantially in their use of the available floral resources, and that this mainly depends on flower attractiveness (floral color and shape). Interestingly, we also detected floral resource partitioning at the colony level within all commonest bumblebee species evaluated. In general, floral resource partitioning between bumblebee individuals decreased with species- and individual-level variation in body size (intertegular span). These results suggest that bumblebee species may coexist via the flexibility in their preferences for specific floral traits, which filters up to support the co-occurrence of high numbers of species and individuals in this global hotspot of species richness.

Designing and Evaluating a Portable UV-LED Vane Trap to Expedite Arthropod Biodiversity Discovery.

A novel design of a portable funnel light trap (PFLT) was presented for collecting insects in ecological studies. The trap consists of a compact plastic box equipped with a light source and power source, along with two plastic polypropylene interception vanes. The PFLT costs 18.3 USD per unit and weighs approximately 300 g. A maximum of six PFLT units can be packed in one medium-sized backpack (32 cm × 45 cm × 15 cm, 20 L), making it easier to set up multiple units in remote areas wherein biodiversity research is needed. The low cost and weight of the trap also allows for large-scale deployment. The design is customizable and can be easily manufactured to fit various research needs. To validate the PFLT's efficacy in collecting insects across different habitat types, a series of field experiments were conducted in South Korea and Laos, where 37 trials were carried out. The PFLT successfully collected 7497 insects without experiencing battery issues or damage by rain or wind. Insect compositions and abundances differed across the three sampled habitat types: forests, grasslands, and watersides. This new FLT trap is an important tool for studying and protecting insect biodiversity, particularly in areas wherein conventional light traps cannot be deployed.

Protein S-palmitoylation enhances profibrotic signaling in response to cadmium.

Cadmium (Cd) is a naturally occurring, toxic environmental metal found in foods. Humans do not have an efficient mechanism for Cd elimination; thus, Cd burden in humans increases with age. Cell and mouse studies show that Cd burden from low environmental levels of exposure impacts lung cell metabolism, proliferation signaling and cell growth as part of disease-promoting profibrotic responses in the lungs. Prior integrative analysis of metabolomics and transcriptomics identified the zDHHC11 transcript as a central functional hub in response to Cd dose. zDHHC11 encodes a protein S-palmitoyltransferase, but no evidence is available for effects of Cd on protein S-palmitoylation. In the present research, we studied palmitoylation changes in response to Cd and found increased protein S-palmitoylation in human lung fibroblasts that was inhibited by 2-bromopalmitate (2-BP), an irreversible palmitoyltransferase inhibitor. Mass spectrometry-based proteomics showed palmitoylation of proteins involved in divalent metal transport and in fibrotic signaling. Mechanistic studies showed that 2-BP inhibited palmitoylation of divalent metal ion transporter ZIP14 and also inhibited cellular Cd uptake. Transcription analyses showed that Cd stimulated transforming growth factor (TGF)-ß1 and ß3 expression within 8 h and lung fibrotic markers α-smooth muscle actin, matrix metalloproteinase-2, and collagen 1α1 gene expression and that these effects were blocked by 2-BP. Because 2-BP also blocked palmitoylation of proteins controlled by TGFß1, these results show that palmitoylation impacts Cd-dependent fibrotic signaling both by enhancing cellular Cd accumulation and by supporting post-translational processing of TGFß1-dependent proteins.

Lessons from assembling UCEs: A comparison of common methods and the case of Clavinomia (Halictidae).

Sequence data assembly is a foundational step in high-throughput sequencing, with untold consequences for downstream analyses. Despite this, few studies have interrogated the many methods for assembling phylogenomic UCE data for their comparative efficacy, or for how outputs may be impacted. We study this by comparing the most commonly used assembly methods for UCEs in the under-studied bee lineage Nomiinae and a representative sampling of relatives. Data for 63 UCE-only and 75 mixed taxa were assembled with five methods, including ABySS, HybPiper, SPAdes, Trinity and Velvet, and then benchmarked for their relative performance in terms of locus capture parameters and phylogenetic reconstruction. Unexpectedly, Trinity and Velvet trailed the other methods in terms of locus capture and DNA matrix density, whereas SPAdes performed favourably in most assessed metrics. In comparison with SPAdes, the guided-assembly approach HybPiper generally recovered the highest quality loci but in lower numbers. Based on our results, we formally move Clavinomia to Dieunomiini and render Epinomia once more a subgenus of Dieunomia. We strongly advise that future studies more closely examine the influence of assembly approach on their results, or, minimally, use better-performing assembly methods such as SPAdes or HybPiper. In this way, we can move forward with phylogenomic studies in a more standardized, comparable manner.

Cryptic diversity begets challenges and opportunities in biodiversity research.

How many species of life are there on Earth? This is a question that we want to know but cannot yet answer. Some scholars speculate that the number of species may reach 2.2 billion when considering cryptic diversity and that each morphology-based insect species may contain an average of 3.1 cryptic species. With nearly two million described species, such high estimates of cryptic diversity would suggest that cryptic species are widespread. The development of molecular species delimitation has led to the discovery of a large number of cryptic species, and cryptic biodiversity has gradually entered our field of vision and attracted more attention. This paper introduces the concept of cryptic species, how they evolve, and methods by which they may be discovered and confirmed, and provides theoretical and methodological guidance for the study of hidden species. A workflow of how to confirm cryptic species is provided. In addition, the importance and reliability of multi-evidence-based integrated taxonomy are reaffirmed as a way to better standardize decision-making processes. Special focus on cryptic diversity and increased funding for taxonomy is needed to ensure that cryptic species in hyperdiverse groups are discoverable and described. An increased focus on cryptic species in the future will naturally arise as more difficult groups are studied, and thereby, we may finally better understand the rules governing the evolution and maintenance of cryptic biodiversity.

A globally synthesised and flagged bee occurrence dataset and cleaning workflow.

Species occurrence data are foundational for research, conservation, and science communication, but the limited availability and accessibility of reliable data represents a major obstacle, particularly for insects, which face mounting pressures. We present BeeBDC, a new R package, and a global bee occurrence dataset to address this issue. We combined >18.3 million bee occurrence records from multiple public repositories (GBIF, SCAN, iDigBio, USGS, ALA) and smaller datasets, then standardised, flagged, deduplicated, and cleaned the data using the reproducible BeeBDC R-workflow. Specifically, we harmonised species names (following established global taxonomy), country names, and collection dates and, we added record-level flags for a series of potential quality issues. These data are provided in two formats, "cleaned" and "flagged-but-uncleaned". The BeeBDC package with online documentation provides end users the ability to modify filtering parameters to address their research questions. By publishing reproducible R workflows and globally cleaned datasets, we can increase the accessibility and reliability of downstream analyses. This workflow can be implemented for other taxa to support research and conservation.

Metabolic reprograming and increased inflammation by cadmium exposure following early-life respiratory syncytial virus infection-the involvement of protein S-palmitoylation.

Early-life respiratory syncytial virus (RSV) infection (eRSV) is one of the leading causes of serious pulmonary disease in children. eRSV is associated with higher risk of developing asthma and compromised lung function later in life. Cadmium (Cd) is a toxic metal, widely present in the environment and in food. We recently showed that eRSV re-programs metabolism and potentiates Cd toxicity in the lung, and our transcriptome-metabolome-wide study showed strong associations between S-palmitoyl transferase expression and Cd-stimulated lung inflammation and fibrosis signaling. Limited information is available on the mechanism by which eRSV re-programs metabolism and potentiates Cd toxicity in the lung. In the current study, we used a mouse model to examine the role of protein S-palmitoylation (Pr-S-Pal) in low dose Cd-elevated lung metabolic disruption and inflammation following eRSV. Mice exposed to eRSV were later treated with Cd (3.3 mg CdCl2/L) in drinking water for 6 weeks (RSV+Cd). The role of Pr-S-Pal was studied using a palmitoyl transferase inhibitor, 2-bromopalmitate (BP, 10 µM). Inflammatory marker analysis showed that cytokines, chemokines and inflammatory cells were highest in the RSV+Cd group, and BP decreased inflammatory markers. Lung metabolomics analysis showed that pathways including phenylalanine, tyrosine and tryptophan, phosphatidylinositol and sphingolipid were altered across treatments. BP antagonized metabolic disruption of sphingolipid and glycosaminoglycan metabolism by RSV+Cd, consistent with BP effect on inflammatory markers. This study shows that Cd exposure following eRSV has a significant impact on subsequent inflammatory response and lung metabolism, which is mediated by Pr-S-Pal, and warrants future research for a therapeutic target.

The evolutionary history of bees in time and space.

Bees are the most significant pollinators of flowering plants. This partnership began ca. 120 million years ago, but the uncertainty of how and when bees spread across the planet has greatly obscured investigations of this key mutualism. We present a novel analysis of bee biogeography using extensive new genomic and fossil data to demonstrate that bees originated in Western Gondwana (Africa and South America). Bees likely originated in the Early Cretaceous, shortly before the breakup of Western Gondwana, and the early evolution of any major bee lineage is associated with either the South American or African land masses. Subsequently, bees colonized northern continents via a complex history of vicariance and dispersal. The notable early absences from large landmasses, particularly in Australia and India, have important implications for understanding the assembly of local floras and diverse modes of pollination. How bees spread around the world from their hypothesized Southern Hemisphere origin parallels the histories of numerous flowering plant clades, providing an essential step to studying the evolution of angiosperm pollination syndromes in space and time.

Launching insectphylo.org; a new hub facilitating construction and use of synthesis molecular phylogenies of insects.

The Holy Grail of an Insect Tree of Life can only be 'discovered' through extensive collaboration among taxon specialists, phylogeneticists and centralized frameworks such as Open Tree of Life, but insufficient effort from stakeholders has so far hampered this promising approach. The resultant unavailability of synthesis phylogenies is an unfortunate situation given the numerous practical usages of phylogenies in the near term and against the backdrop of the ongoing biodiversity crisis. To resolve this issue, we establish a new online hub that centralizes the collation of relevant phylogenetic data and provides the resultant synthesis molecular phylogenies. This is achieved through key developments in a proposed pipeline for the construction of a species-level insect phylogeny. The functionality of the framework is demonstrated through the construction of a highly supported, species-comprehensive phylogeny of Diptera, built from integrated omics data, COI DNA barcodes, and a compiled database of over 100 standardized, published Diptera phylogenies. Machine-readable forms of the phylogeny (and subsets thereof) are publicly available at insectphylo.org, a new public repository for species-comprehensive phylogenies for biological research.

The impact of heavy metal pollution on wild bee communities in smallholder farmlands.

Wild bees provide important pollination services, but they face numerous stressors that threaten them and their ecosystem services. Wild bees can be exposed to heavy metal pollution through the consumption of nectar, pollen, and water, which might cause bee decline. While some studies have measured heavy metal concentrations in honeybees, few studies have monitored heavy metal concentrations in wild bees or explored their potential effects on wild bee communities. To investigate the impact of heavy metal pollution on wild bee communities, heavy metal concentrations, including vanadium (V), chromium (Cr), nickel (Ni), cadmium (Cd), Zinc (Zn) and lead (Pb) in multiple wild bee species were measured. Multiple wild bee species, including: Xylocopa tranquabaroroum, Eucera floralia, Apis cerana, and small bee mixtures (representing multiple small wild bee species) were sampled from 18 sites in Quzhou, Zhejiang Province, China. The findings demonstrated that there were significant differences in heavy metal concentrations among different bee species. The concentrations of V, Zn, Cd, and Pb in X. tranquabaroroum, the largest bee species in this study, were lower than that in the other three sample groups. Furthermore, there were significant negative correlations between heavy metal pollution and wild bee diversity and species richness, but not with abundance. Particularly, there was no significant relationship between heavy metal pollution and the abundance of small bees. Given these worrying findings, monitoring multiple heavy metals in wild bees should be conducted for protecting wild bee diversity and securing their pollination services.

Substrate-dependent metabolomic signatures of myeloperoxidase activity in airway epithelial cells: Implications for early cystic fibrosis lung disease.

Myeloperoxidase (MPO) is released by neutrophils in inflamed tissues. MPO oxidizes chloride, bromide, and thiocyanate to produce hypochlorous acid (HOCl), hypobromous acid (HOBr), and hypothiocyanous acid (HOSCN), respectively. These oxidants are toxic to pathogens, but may also react with host cells to elicit biological activity and potential toxicity. In cystic fibrosis (CF) and related diseases, increased neutrophil inflammation leads to increased airway MPO and airway epithelial cell (AEC) exposure to its oxidants. In this study, we investigated how equal dose-rate exposures of MPO-derived oxidants differentially impact the metabolome of human AECs (BEAS-2B cells). We utilized enzymatic oxidant production with rate-limiting glucose oxidase (GOX) coupled to MPO, and chloride, bromide (Br-), or thiocyanate (SCN-) as substrates. AECs exposed to GOX/MPO/SCN- (favoring HOSCN) were viable after 24 h, while exposure to GOX/MPO (favoring HOCl) or GOX/MPO/Br- (favoring HOBr) developed cytotoxicity after 6 h. Cell glutathione and peroxiredoxin-3 oxidation were insufficient to explain these differences. However, untargeted metabolomics revealed GOX/MPO and GOX/MPO/Br- diverged significantly from GOX/MPO/SCN- for dozens of metabolites. We noted methionine sulfoxide and dehydromethionine were significantly increased in GOX/MPO- or GOX/MPO/Br--treated cells, and analyzed them as potential biomarkers of lung damage in bronchoalveolar lavage fluid from 5-year-olds with CF (n = 27). Both metabolites were associated with increasing bronchiectasis, neutrophils, and MPO activity. This suggests MPO production of HOCl and/or HOBr may contribute to inflammatory lung damage in early CF. In summary, our in vitro model enabled unbiased identification of exposure-specific metabolite products which may serve as biomarkers of lung damage in vivo. Continued research with this exposure model may yield additional oxidant-specific biomarkers and reveal explicit mechanisms of oxidant byproduct formation and cellular redox signaling.

Low-dose vanadium pentoxide perturbed lung metabolism associated with inflammation and fibrosis signaling in male animal and in vitro models.

Vanadium is available as a dietary supplement and also is known to be toxic if inhaled, yet little information is available concerning the effects of vanadium on mammalian metabolism when concentrations found in food and water. Vanadium pentoxide (V+5) is representative of the most common dietary and environmental exposures, and prior research shows that low-dose V+5 exposure causes oxidative stress measured by glutathione oxidation and protein S-glutathionylation. We examined the metabolic impact of V+5 at relevant dietary and environmental doses (0.01, 0.1, and 1 ppm for 24 h) in human lung fibroblasts (HLFs) and male C57BL/6J mice (0.02, 0.2, and 2 ppm in drinking water for 7 mo). Untargeted metabolomics using liquid chromatography-high-resolution mass spectrometry (LC-HRMS) showed that V+5 induced significant metabolic perturbations in both HLF cells and mouse lungs. We noted 30% of the significantly altered pathways in HLF cells, including pyrimidines and aminosugars, fatty acids, mitochondrial and redox pathways, showed similar dose-dependent patterns in mouse lung tissues. Alterations in lipid metabolism included leukotrienes and prostaglandins involved in inflammatory signaling, which have been associated with the pathogenesis of idiopathic pulmonary fibrosis (IPF) and other disease processes. Elevated hydroxyproline levels and excessive collagen deposition were also present in lungs from V+5-treated mice. Taken together, these results show that oxidative stress from environmental V+5, ingested at low levels, could alter metabolism to contribute to common human lung diseases.NEW & NOTEWORTHY We used relevant dietary and environmental doses of Vanadium pentoxide (V+5) to examine its metabolic impact in vitro and in vivo. Using liquid chromatography-high-resolution mass spectrometry (LC-HRMS), we found significant metabolic perturbations, with similar dose-dependent patterns observed in human lung fibroblasts and male mouse lungs. Alterations in lipid metabolism included inflammatory signaling, elevated hydroxyproline levels, and excessive collagen deposition were present in V+5-treated lungs. Our findings suggest that low levels of V+5 could trigger pulmonary fibrotic signaling.

Phylogeny-based assignment of functional traits to DNA barcodes outperforms distance-based, in a comparison of approaches.

The full potential for using DNA barcodes for profiling functional trait diversity has yet to be determined in plants and animals; thus, we outline a general framework for quantifying functional trait diversity of insect community DNA and propose and assess the accuracy of three methods for achieving this. We built a novel dataset of traits and DNA barcodes for wild bees in China. An informatics framework was developed for phylogeny-based integration of these data and prediction of traits for any subject barcodes, which was compared with two distance-based methods. For Phylogenetic Assignment, we additionally conducted a species-level analysis of publically available bee trait data. Under the specimen-level dataset, the rate of trait assignment was negatively correlated with distance between the query and the nearest trait-known reference, for all methods. Phylogenetic Assignment was found to perform best under several criteria; particularly, it had the lowest false-positive rate (rarely returning a state prediction where success was unlikely; where the distance from query to the nearest reference was high). For a wider range of compiled traits, conservative life-history traits showed the highest rates of assignment; for example, sociality was predicted with confidence at 53%, parasitism at 44% and nest location at 33%. As outlined herein, automated trait assignment might be applied at scale to either barcodes or metabarcodes. With further compilation and databasing of DNA barcode and trait data, the rate and accuracy of trait assignment is expected to increase to the point of being a widely viable and informative approach.

Renin-Angiotensin System Modulation With Synthetic Angiotensin (1-7) and Angiotensin II Type 1 Receptor-Biased Ligand in Adults With COVID-19: Two Randomized Clinical Trials.

Importance: Preclinical models suggest dysregulation of the renin-angiotensin system (RAS) caused by SARS-CoV-2 infection may increase the relative activity of angiotensin II compared with angiotensin (1-7) and may be an important contributor to COVID-19 pathophysiology. Objective: To evaluate the efficacy and safety of RAS modulation using 2 investigational RAS agents, TXA-127 (synthetic angiotensin [1-7]) and TRV-027 (an angiotensin II type 1 receptor-biased ligand), that are hypothesized to potentiate the action of angiotensin (1-7) and mitigate the action of the angiotensin II. Design, Setting, and Participants: Two randomized clinical trials including adults hospitalized with acute COVID-19 and new-onset hypoxemia were conducted at 35 sites in the US between July 22, 2021, and April 20, 2022; last follow-up visit: July 26, 2022. Interventions: A 0.5-mg/kg intravenous infusion of TXA-127 once daily for 5 days or placebo. A 12-mg/h continuous intravenous infusion of TRV-027 for 5 days or placebo. Main Outcomes and Measures: The primary outcome was oxygen-free days, an ordinal outcome that classifies a patient's status at day 28 based on mortality and duration of supplemental oxygen use; an adjusted odds ratio (OR) greater than 1.0 indicated superiority of the RAS agent vs placebo. A key secondary outcome was 28-day all-cause mortality. Safety outcomes included allergic reaction, new kidney replacement therapy, and hypotension. Results: Both trials met prespecified early stopping criteria for a low probability of efficacy. Of 343 patients in the TXA-127 trial (226 [65.9%] aged 31-64 years, 200 [58.3%] men, 225 [65.6%] White, and 274 [79.9%] not Hispanic), 170 received TXA-127 and 173 received placebo. Of 290 patients in the TRV-027 trial (199 [68.6%] aged 31-64 years, 168 [57.9%] men, 195 [67.2%] White, and 225 [77.6%] not Hispanic), 145 received TRV-027 and 145 received placebo. Compared with placebo, both TXA-127 (unadjusted mean difference, -2.3 [95% CrI, -4.8 to 0.2]; adjusted OR, 0.88 [95% CrI, 0.59 to 1.30]) and TRV-027 (unadjusted mean difference, -2.4 [95% CrI, -5.1 to 0.3]; adjusted OR, 0.74 [95% CrI, 0.48 to 1.13]) resulted in no difference in oxygen-free days. In the TXA-127 trial, 28-day all-cause mortality occurred in 22 of 163 patients (13.5%) in the TXA-127 group vs 22 of 166 patients (13.3%) in the placebo group (adjusted OR, 0.83 [95% CrI, 0.41 to 1.66]). In the TRV-027 trial, 28-day all-cause mortality occurred in 29 of 141 patients (20.6%) in the TRV-027 group vs 18 of 140 patients (12.9%) in the placebo group (adjusted OR, 1.52 [95% CrI, 0.75 to 3.08]). The frequency of the safety outcomes was similar with either TXA-127 or TRV-027 vs placebo. Conclusions and Relevance: In adults with severe COVID-19, RAS modulation (TXA-127 or TRV-027) did not improve oxygen-free days vs placebo. These results do not support the hypotheses that pharmacological interventions that selectively block the angiotensin II type 1 receptor or increase angiotensin (1-7) improve outcomes for patients with severe COVID-19. Trial Registration: ClinicalTrials.gov Identifier: NCT04924660.

Antagonistic Interactions in Mitochondria ROS Signaling Responses to Manganese.

Antagonistic interaction refers to opposing beneficial and adverse signaling by a single agent. Understanding opposing signaling is important because pathologic outcomes can result from adverse causative agents or the failure of beneficial mechanisms. To test for opposing responses at a systems level, we used a transcriptome-metabolome-wide association study (TMWAS) with the rationale that metabolite changes provide a phenotypic readout of gene expression, and gene expression provides a phenotypic readout of signaling metabolites. We incorporated measures of mitochondrial oxidative stress (mtOx) and oxygen consumption rate (mtOCR) with TMWAS of cells with varied manganese (Mn) concentration and found that adverse neuroinflammatory signaling and fatty acid metabolism were connected to mtOx, while beneficial ion transport and neurotransmitter metabolism were connected to mtOCR. Each community contained opposing transcriptome-metabolome interactions, which were linked to biologic functions. The results show that antagonistic interaction is a generalized cell systems response to mitochondrial ROS signaling.

The Whole Genome of Colletes collaris (Hymenoptera: Colletidae): An Important Step in Comparative Genomics of Cellophane Bees.

Gaining knowledge on bees is of the utmost importance due to the paramount role that they play in angiosperm pollination. Herein, we provide the first genome assembly of Colletes collaris, a pan-Eurasian cellophane bee. We sequenced 50.53 Gbp of long-read data plus 57.36 Gbp of short-read data in Oxford Nanopore Technologies and Illumina platforms, respectively. The genome assembly consisted of 374.75 Mbp distributed across 374 contigs, with L50 and N50 of 9 and 8.96 Mbp, respectively. We predicted the genome to comprise 20,399 protein-coding genes, 467,947 repeats, and 4,315 non-coding RNA genes. The transcriptome and mitochondrial genome of the species were also assembled. Gene family analysis with 15 insect species identified 14,417 families, 9,517 of them found in C. collaris. A dated phylogenomic analysis revealed high numbers of orthogroups experiencing rapid evolution within Colletes.